function a=rephase(b)
c=length(b);a=zeros(1,c);a(1)=b(1);e=0;
for d=2:c
    if  b(d)-b(d-1)>pi
        e=e-2*pi;
    elseif b(d)-b(d-1)<-pi
        e=e+2*pi;
    end
    a(d)=b(d)+e;
end